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Abstract. We present the results of a BeppoSAX ob- 
servation of PSR J0218-h4232 which, for the first time, 
provides detailed information on the pulsar's temporal 
and spectral emission properties in the broad band 1- 
10 keV. We detected X-ray pulses with a pulsed frac- 
tion of 73±12 %. The pulse profile is characterized by 
two peaks phase separated by Acj) = 0.47 ± 0.05. The 
pulsed spectrum is best described by a power-law of pho- 
ton index 0.61±0.32 with an unabsorbed (2-10 keV) X-ray 
flux of 4.1x10"-'^^ erg cm~^ s~^ implying a luminosity of 
= 1.3 X 10^2 e (d/5.7 kpc)2 erg s'^ and an X-ray ef- 
ficiency of L^/E = 4.8 X 10-4 e (d/5.7 kpc)^ where 9 is 
the solid angle spanned by the emission beam. 
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Stars: neutron. Pulsar: individual: PSR J0218+4232, X- 
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1. Introduction 

PSR J0218+4232 is a 2.3 ms pulsar in a two day orbit 
around a ~ 0.2 Mq white dwarf companion (Navarro et 
al. 1995) with a period derivative of P = 8 x 10~^° s s~^. 
The pulsar has a spin-down energy of £' = 2.5 x 10^^ erg 
s~^, a dipolar magnetic field component at the star surface 
of B± = 4.3 x 10^ G and a spin-down age of < 4.6 x 
10® years. The pulsar distance inferred from its dispersion 
measure and from the electron density model of Taylor & 
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Cordes (1993) is > 5.7 kpc. A detection of the companion 
star at optical wavelength was reported recently by Van 
Kerkwijk (1996). 

Soft X-ray emission from the pulsar was first detected 
by Verbunt et al. (1996) using the ROSAT HRI. A follow- 
up observation confirmed the detection and discovered X- 
ray pulsation at a significance of about 5 cr using 200 HRI 
counts (Kuiper et al. 1998). The X-ray pulse profile de- 
duced from those data is characterized by a sharp main 
pulse with an indication of a second peak at a phase sep- 
aration of A(f> 0.47. The pulsed fraction inferred from 
the ROSAT HRI data is 37±13 %. Furthermore, Kuiper 
et al. (1998) show that the measured large DC compo- 
nent is consistent with a ~ 14" diameter compact nebula 
surrounding the pulsar, but confirmation is required. How- 
ever, it is interesting to note that also in the radio domain 
the source exhibits an unusually high DC component of 
^ 50 % (Navarro et al. 1995). The HRI provides no spec- 
tral information and the number of counts recorded in a 
serendipitous off-axis PSPC observation does not allow 
spectral modeling. Also ASCA detected this source, how- 
ever, the observation was too short: no pulsation could be 
seen, and a spectral fit with a power-law photon index 
of 1.6±0.6 could only be made to the weak total excess 
(Kawai & Saito 1999). Therefore, no detailed spectral in- 
formation on the pulsed X-ray emission was available prior 
to the BeppoSAX observation reported in this paper (see 
Becker & Triimper 1999 for a review of the X-ray proper- 
ties of millisecond pulsars). 

Noting the spatial coincidence of PSR J0218+4232 
with the EGRET source 2EG J0220+4228, Verbunt et 
al. (1996) tentatively identified the pulsar with the high- 
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energy 7-ray source. Using some additional EGRET ob- 
servations, and applying a combination of spatial and 
timing analysis, Kuiper et al. (1999) conclude that 2EG 
J0220+4228 is probably a multiple source: between 0.1 
and 1 GeV PSR J0218-I-4232 is the most likely counter- 
part, and above 1 GeV the bright BL Lac 3C 66 A is the 
best candidate counterpart. The third EGRET catalog 
(Hartman et al. 1999), which is based on more viewing 
periods than the 2EG catalog, identifies 3EG J0222-H4253 
(2EG J0220+4228) with 3C 66A, rather than with the ms- 
pulsar. However, in a note on this source, they indicate 
that the identification with 3C 66A stems from the cat- 
alog position based on the > 1 GeV map. Furthermore, 
they confirm that for lower energies (100-300 MeV) the 
EGRET map is consistent with all the source flux coming 
from the pulsar, 3C 66A being statistically excluded. 

In this paper we present the results of a BeppoSAX 
observation of PSR J0218-f 4232 which, for the first time, 
provides detailed information on the pulsar's temporal and 
spectral emission properties across a wide band from 1 to 
10 keV. 

2. Observation 

PSR J0218-h4232 was observed on January 1999 14-16 by 
the Narrow Field Instruments (NFIs) aboard BeppoSAX 
satellite (Boclla et al. 1997a). We report in this paper 
only results from data collected with the MEGS instru- 
ment (Boella et al. 1997b), sensitive in the energy range 
1-10 keV; the net exposure time was 82 795 s. The LEGS 
instrument (0.1-10 keV; Parmar et al. 1997) observed the 
source for a much shorter time due to the constraint on 
operating only during spacecraft night. The HPGSPG (4- 
60 keV; Manzo et al. 1997) was not working during the 
observation. A signal detected with the PDS (13-200 keV; 
Frontera et al. 1997) will be shortly discussed in Section 
5. No pulsed signals have been detected in the LEGS and 
PDS data. 

During the observation only two of the three MEGS de- 
tectors were operating (MEGS2 and MEGS3). MEGS has 
a field of view of 56' and an angular resolution of about 
1.'2 at 6 keV. During the observation of PSR J0218+4232 
the instrument operated in direct mode, transmitting to 
ground information on each individual photon. Standard 
procedures and selection criteria were applied on the ob- 
servation data to avoid the South Atlantic Anomaly, solar, 
bright Earth and particle contaminationaH. Event reduc- 
tion has been performed using the SAXDAS v. 2. 0.0 pack- 
age. 

3. Spatial analysis 

A 100 ks ROSAT HRI (0.1-2.4 keV) observation of PSR 
J0218+4232 (see Kuiper et al. 1998) revealed 7 X-ray 



sources within a radius of ^ 5' around our target (see Fig. 
1). Given the extended tails of the MEGS Point-Spread 
Function (PSF), accurate spatial analysis is required to 
separate the signal of PSR J0218-f 4232 from that of any 
neighbouring sources. We used a Maximum Likelihood ap- 
proach in searching for individual sources on top of a back- 
ground model, as well as for analyzing simultaneously sev- 
eral sources and a background model. 
The search is applied on various equidistant trial posi- 
tions within a selected part of the instrument field of view, 
where it checks for the presence of a source taking into ac- 
count the Poissonian nature of the data. Maximum Like- 
lihood Ratio (MLR) tests are performed describing in the 
zero Ti.Q hypothesis the 2-d event distribution in terms of 
a fiat background model only, while in the alternative 7ii 
hypothesis the data is described in terms of a point source 
at the trial position and a fiat background model (see e.g. 
Kuiper et al. 1998 for a more detailed description). 

The improvement in the likelihood C between Tio and 
Hi computed from the quantity Q = -2\n{C"° /C^^) 
yields simultaneously the detection significance and 
strength of the source. The probability distribution of Q 
is that of a for ni — uq degree of freedom (d.o.f.) where 
riQ and rii are the d.o.f. for the Tio a-nd Tii hypotheses, re- 
spectively. In searching for sources ni — uq — 3, while in a 
detection test of a known source ni—riQ = 1. The distribu- 
tion of Q as a function of trial position is the MLR-map. 

We applied the above approach to the MEGS data 
without energy selection allowing for maximum source 
statistics. We used an energy averaged PSF assuming a 
power-law spectral shape of index 1.5, given by the fol- 
lowing expression: 
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PSF{x, y) = 
' E-^-^ ■ {S2{E) + Sz{E)) ■ PSF{x,v I E)-dE 



tittp:/ /www. sdc.asi.it/software/cookbook as a reference 



about the data analysis software and reduction procedures. 



E~^-^ ■ {S2{E) + S:,{E)) ■ dE ' 

where S2{E) and Ss{E) are the sensitive areas of the two 
operating MEGS units. The PSF{x,y | E) is modelled 
using the parametric form given in Boella et al. (1997b). 
Furthermore, we investigated the impact of the assumed 
power-law index on the MLR results and found that this 
is negligible, as expected by the moderate dependence on 
energy of the PSF. 

The resulting MLR map/image is shown in Fig. la. 
The maximum value in this figure (> 21 ct) is reached at a 
position consistent with that determined by the ROSAT 
HRI for PSR J0218-I-4232 (indicated by a star symbol). 
It is, however, evident that several of the sources detected 
by ROSAT below 2.4 keV (indicated by filled triangles) 
contribute significantly to the extended excess. We can 
show this better, by repeating the search for sources on 
top of a fiat background model and a point source at 
the pulsar position (effectively "subtracting" the counts 
from PSR J0218-h4232 from Fig. la). The remaining ex- 
tended feature (Fig. lb) clearly follows the distribution of 
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Fig. 1. Maximum Likelihood Ratio maps in instrumental coordinates using all MECS data (~ 1-10 keV) for the 

sky region around PSR J0218+4232 (star). The contours start at a value of 3 ct source detection significance with a stepsize 
of 2 0" (for 1 d.o.f.) or equivalently at Q values (see text) of 9, 25, 49, . . .. The positions of all sources detected by the ROSAT 
HRI observation below 2.4 keV are indicated by filled triangles, a) Evidence for sources above a flat background showing PSR 
J0218+4232 (~ 21 a) on top of an extended excess, b) Like a), but PSR J0218+4232 has now been "subtracted", c) No significant 
excesses are left when 9 ROSAT sources are added to the background model 



the nearby sources, reaching detection significances up to 
about 10 cr. Figure Ic shows the MLR map resulting from 
a point source search on top of a flat background and 9 
sources fixed at the ROSAT HRI source positions: no sig- 
nificant residual emission remains. The total number of 
counts assigned by this analysis to PSR J0218+4232 is 
403±29. 

The position of PSR J0218+4232 determined from the 
MLR map shown in Fig. la is shifted by ~ 23" with re- 
spect to the radio position, well within the systematic un- 
certainty of the accuracy of the BeppoSAX pointing re- 
construction H. 

The spatial analysis allows also the evaluation of an op- 
timum event extraction radius by computing the Signal- 
to-Noise ratio as a function of radial distance r from the 
source position. Taking as source radial profile the PSF, 
and as noise the measured radial count distribution, we ob- 
tained the value r — 100". Note that none of the nearby 
sources fall within this radius. The selected region con- 
tains 62 % (250 counts) of the MECS source signal; the 
number of background events in the same area is 92 with 
a contribution from the neighbouring sources of 12 counts. 

4. Timing analysis 

The arrival times of all selected events were converted to 
the Solar System Barycentric Frame using BARYCONVQ, 
then folded according to the radio ephemeris (see Table 1), 
correcting for the pulsar binary motion. The pulse phase 
distribution deviates from a statistically flat distribution 



at a 6.8 a level applying a test (Buccheri et al. 1983) 
using the first two harmonics. 

The MECS timing resolution is dominated by the pre- 
cision of the OBT (On Board Time)-UTC conversion: an 
inspection of the residuals of the OBT-UTC linear fit re- 
veals a systematic scatter with an rms of the order of ~ 0.2 
ms that cannot be reduced by fitting higher order poly- 
nomials. This effect corresponds to an uncertainty in the 
photon phases of 0.087. 



Table 1. Ephemeris of PSR J0218-H4232 



Parameter 



Value 



Right Ascension (J2000) 02*^ 18™ 6^:350 

Declination (J2000) 42° 32' 17'.'44 

Epoch validity start/end (MJD) 49092 - 50900 

Frequency 430.4610674213 Hz 



Frequency derivative 

Epoch of the period (MJD) 

Orbital period 

apsin i 

Eccentricity 

Longitude of periastron 

Time of ascending node (MJD) 



-1.4342 X 10~'* Hz s" 

49996.000000023 

175292.3020 s 

1.98444 (h-s) 





49996.637640 



http:/ /www. sdc.asi.it/software/cookbook/attitude.htmj 

see pittp:/ /www. sdc.asi.it/software/saxdas/baryconv. html 



The light curve resulting from folding all MECS events 
(1-10 keV) with a phase resolution of 18 bins 0.13 
ms) is shown in Fig. 2 (top panel): two cycles are plotted 
for clarity. The MECS background level, indicated with 
a grey band (±1 cr errors), is determined in the spatial 
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Fig. 2. PSR J0218+4232 light curve in the 1-10 keV range as 
observed by the MECS (top panel) and at 610 MHz (bottom 
panel, Stairs et al. 1999). In the top panel a typical statistical 
1 cr error bar is shown and the MECS background level as 
determined in the spatial analysis is indicated as a grey band, 
representing the ±1 a widths. The two horizontal bars above 
the MECS profile indicate the definitions of the phase intervals 
for the two peaks. Peak 1 (PI) and Peak 2 (P2), as explained 
in the text. The alignment between the two curves is arbitrary 



analysis and therefore includes the contributions from the 
neighbouring sources. In the same figure (bottom panel) 
the 610 MHz profile as derived by Stairs et al. (1999) is 
also presented. The phase relation between the two pro- 
files is unknown for the time being: they have been arbi- 
trarily aligned. Notice that, like in the radio light curve, 
the MECS profile is not simple: it has a prominent dou- 
ble peak structure with a relative phase separation (centre 
to centre) of 0.47±0.05, confirming the value derived by 
Kuiper et al. (1998) from ROSAT data below 2.4 keV. 
The phase histograms for 1.6-4 keV and 4-10 keV are 



shown respectively in the middle and bottom panels of 
Fig. 3. The ROSAT (0.1-2.4 keV) profile (Kuiper et al. 
1998), shifted in phase to obtain the highest peak coin- 
cident with the most significant one in the MECS softer 
light curve (middle panel), is also shown for comparison 
in the top panel of Fig. 3. These profiles clearly show a 
change of the relative peak intensities. The peak at phase 
0.8 is stronger in the low energy histograms (top and mid- 
dle panels), while that at phase 0.3 is more prominent at 
high energies (bottom panel). 

The structure measured in the MECS light curves of 
two peaks separated by narrow valleys and the uncertainty 
due to the OBT-UTC conversion residuals make it diSi- 
cult to establish an unpulsed interval: slightly different 
values could affect the unpulsed level. Notice, moreover, 
that while in the ROSAT profile a DC component is appar- 
ent above the background level, the same does not hold 
for the 4-10 keV profile, where the background level is 
consistent with those measured in the valleys of the light 
curve. In the intermediate energy profile, at face value, 
there might be some evidence of a DC component but not 
as high as in the ROSAT profile. Such effect could be due 
to the presence of a quite soft unpulsed emission. In an 
attempt to quantify the evidence for a DC component, 
we applied the bootstrap method proposed by Swanepoel 
et al. (1996), which allows to estimate the DC level, and 
therefore also the pulsed fraction, working directly on the 
phases of individual photons. However, in the available 
form, this method is able to find only one unpulsed inter- 
val and cannot account for systematic errors. One should 
therefore realize that the quoted errors on the calculated 
parameters are only statistical. Using this algorithm we 
found, for the whole MECS energy band, an unpulsed in- 
terval of 0.47-0.61 and a corresponding pulsed fraction of 
0.73±0.12 (1 a error). Notice that the average count level 
in the interval —0.08-0.15 is statistically not different from 
that measured in the former interval. The resulting DC- 
fractions for the three bands of Fig. 3 are shown in Fig. 4. 
The decrease with energy is clearly visible, but the error 
bars are large and the trend deviates only at the 2.4 u 
level from a flat distribution. Given the large error bars 
and the systematic uncertainties, we did not derive a de- 
tailed spectrum of the DC component. 

The relatively high scatter in the residuals of our tim- 
ing solution (up to ~0.2 ms) does not allow to determine 
accurately the intrinsic widths of the two pulses. The total 
scatter in the event arrival times at the SSB might even 
be consistent with the measured width of PI in Figure 2. 

We defined the phase intervals containing the two 
pulses as follows: Peak 1 (PI) phase interval 0.61-0.95 
and Peak 2 (P2) 0.10-0.47, where PI is identified with the 
obvious ROSAT peak and appears to have a softer spec- 
trum than P2. The statistical significance of the suggested 
change in morphology of the phase histogram with energy 
can be evaluated by the P2/P1 ratio. This ratio, computed 
in the energy bands 1.6-4 keV and 4-10 keV after back- 



T. Mineo et al.: Pulse shape and spectrum of PSR J0218+4232 observed with BeppoSAX 



5 



30 



20 



10 




27 



18 



30 



20 



10 



0.1-2.4 keV 



1.6-4 keV 




4-10 keV 



J 



n 



n 



0.0 0.5 1.0 1.5 

Pulse phase 



2.0 



Fig. 3. PSR J0218+4232 phase histograms: energy range 
0.1-2.4 keV from the ROSAT HRI (Kuiper et al. 1998); energy 

ranges 1.6-4 koV and 4-10 kcV from the BeppoSAX MEGS 
(this work). Note the increase of the P2 peak strength with 
respect to PI with energy. The hatched area represents the 
unpulsed level (±1 a) determined in a bootstraping analysis of 
the phase histogram, while the grey area represents the back- 
ground level with 1 a errors determined in the spatial analysis 



ground subtraction, is shown in Fig. 5. The ROSAT point 

derived from the data of Kuiper et al. (1998), indicated 
by a filled circle, is evaluated using the same phase defini- 
tions. Including the ROSAT point, a correlation of P2/P1 
with energy is evident, but the statistical significance is 
only 2 cr, and confirmation is needed. 
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Fig. 4. The DC-fraction as a function of energy. The ROSAT 
point (filled circle) was derived using the same width of the 
unpulsed interval as obtained for the whole MECS range 



5. Spectral analysis 

Spectral fits were performed on the background subtracted 
MECS data binned logarithmically after checking that 
each single energy channel contains at least 20 photons. 
The Nh was fixed to the value of 5x10^° cm~^ (see Ver- 
bunt et al. 1996). We fitted the total pulsed spectrum, i.e. 
the excess counts above the unpulsed level in the phase 
histogram, with a power-law model, taking into account 
the energy dependent interstellar absorption. The result- 
ing photon index is a = 0.61±0.32, with a 2-10 keV fiux 
F = 4.1 X 10~^^ erg cm~^ s~^, and a reduced Xr of 
0.7 (3 dof). A fit with a black body spectral distribution 
gives kT = 2.3 ± 0.5 keV, with a Xr = 1-4 (3 dof). Both 
the values are acceptable, the black body temperature 
T = 2.6 X 10^ K, however, is quite higher than the val- 
ues measured for the thermal components of other pulsars 
that typically does not exceed ^ 5 x 10^ K (Greiveldinger 
et al. 1996, Becker & Triimpcr 1997, 1999). Fig. 6 also 
shows the spectra for the pulses PI and P2 separately, us- 
ing the same energy bins as for the total pulsed spectrum. 
We fitted again power-law models for each of the spectra. 
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but the low statistics result in large errors in the estimated 

parameters. The measured spectral indices of PI and P2 
are 0.84 ± 0.35 and 0.42 ± 0.36 respectively, in agreement 
with the trend seen in the P2/P1 ratio. 

Finally, we determined the total spectrum (pulsed plus 
DC) applying the Maximum Likelihood approach to the 
spatial distributions in smaller energy intervals (see also 
Fig. 6). The fit with a power-law gives a photon spectral 
index a = 0.94±0.22 with a 2-10 keV flux F = 4.3 x lO'^^ 
erg cm~^ s~^, about 10 % higher than the pulsed value 
given above, and a reduced Xr of 0.56 (8 dof). 
Given the hard spectra found above, we also analysed the 
PDS data to investigate whether there was any signal from 
PSR J0218+4232 at higher energies, even though an ex- 
trapolation of the total spectrum measured in the MECS 
to the PDS range predicts a flux well below the PDS sen- 
sitivity threshold. Indeed, no pulsed signal was found, but 
in the 17-25 keV band a DC signal was detected at the 
~4 fj level. Moreover the fit of MECS+PDS spectra with 
a power-law plus a constant factor to take into account 
the intercalibration uncertainties between the two instru- 
ments leads to a spectral index of 1.4 and to a value for 
the constant ^8, well above the expected range of variabil- 
ity (0.7-1, Cusumano et al. 2000). Furthermore, the PDS 
light curve shows that most of the source counts are con- 
centrated in the first half of the observation. Therefore, 
it is likely that we observed a variable source in the large 
field of view of the PDS coUimators (1?3 x 1?3 FWHM). 
Note that there are several sources within a 3° region in 
the soft X-ray catalogs, but none in the hard X-ray ones. 

6. Discussion 

The BeppoSAX observation of PSR J0218-I-4232 provided 

the first detection of pulsed emission from this millisecond 
pulsar for energies up to 10 keV. The source shows a dou- 
ble peaked pulse shape with a remarkably flat (photon 
index ^0.6) spectral distribution and indications that the 
peak intensity ratio depends on photon energy. In partic- 
ular, the P2 peak, which is not prominent in the ROSAT 
low energy phase histogram, becomes the dominant fea- 
ture above 5 keV. 

Our fit for a black body spectrum gives the quite high 
value of 2.6 x 10^ K for the temperature of the emitting 
region. If we assume that the radiation comes from the 
pulsar polar cap, heated, for instance, by the interactions 
with starward moving high energy particles, we can com- 
pute the expected luminosity. The polar cap area of PSR 
J0218-I-4232, defined as usual by the open field lines, is of 
the order of ^ ~ 27ri?2[l - ^(1 - RQ/c)] = 2.9 x 10" 
cm^, where R is the stellar radius; from the Stefan's law 
we obtain a bolometric luminosity of 5 x 10^^ erg s~^, a 
value that must be doubled if the other cap is taken into 
account. Such luminosity is several orders of magnitude 
higher than that derived from the phase averaged flux 
Lx = 1.3 X 10^^ 6 {d/5.7 kpc)^ erg s~^ (corresponding to 



an X-ray efficiency of L^/E = 4.8 x 10"* 6 {d/5.7 kpc)^) 

where Q is the solid angle spanned by the emission beam. 
It is moreover significantly greater than the spin-down en- 
ergy loss rate. This inconsistency could be solved with the 
assumption that the emission spot covers a very small frac- 
tion, of the order of about lO^'' or even less, of the whole 
cap. The flux of heating particles should be then colli- 
mated within a very narrow angle of 0^.2 or smaller. 
Furthermore, the heating flux into the two polar caps must 
be different to explain the spectral difference of the two 
peaks. A non-thermal origin of the X-ray emission seems 
therefore more convincing. 

The variation of the peak ratio with energy recalls a 
similar effect observed for the Crab pulsar (see, for in- 
stance, the recent BeppoSAX data described by Mineo et 
al. 1997). The possibility of detecting such effect could 
be due to a particular orientation among the rotation 
and magnetic axes and the line of sight. The Crab has 
an orthogonal alignment, however, Navarro et al. (1995), 
on the basis of the radio pulse shape, inferred for PSR 
J0218-I-4232 a nearly aligned rotator. The same conclu- 
sion is reached by Stairs et al. (1999) who analyzed high- 
precision, coherently-dispersed polarization profiles at the 
two frequencies of 410 and 610 MHz. Their Rotating Vec- 
tor Model (RVM) fits support the classification of PSR 
J0218+4232 as a nearly aligned rotator with magnetic in- 
clination consistent with O'', namely (8±11) deg at 410 
MHz and (8±15) at 610 MHz. Unfortunately, the impact 
angle parameters for their RVM fits have large uncertain- 
ties and therefore the line-of-sight inclination is uncon- 
strained. 
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Fig. 5. The energy dependence of P2/P1, according to the 
phase interval defined in the text. The ROSAT point is indi- 
cated by a filled circle 
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The Crab-like double pulse X-ray profile of PSR 

J02 18+4232 and its hard non-thcrnial spectrum suggests 
a common magnetospheric origin of the high-energy emis- 
sion. The comparison with Crab is additionally enforced 
by the similar values of the magnetic field at the Ught 
cylinder, as already pointed out by Kuiper et al. (1998) 
and Kawai & Saito (1999). This follows directly from the 
classical assumption that the pulsar braking is entirely 
due to the Poynting flux at the light cylinder, indepen- 
dently of the actual structure of the magnetic field in the 
magnetosphere: 



(2) 



where / is the moment of inertia. For the two pulsars we 
obtain then: 

BL = Bj^M-^)H^)-H^f\ (3) 

J^Cr ^Cr Pcr 

where the suffix Cr is referred to Crab. Assuming that the 
moments of inertia of the two pulsars are equal, we found 
Bl ~ \BL,Cr ^ 3 X 10^ Gauss. 

The nature of the non-thermal mechanism responsible 
for the X-ray emission is still unclear. A possibility is that 
X-ray photons are produced by curvature radiation from 
relativistic electrons moving outwards along the field lines 
in proximity of the light cylinder. The Lorentz factor of 
these particles to radiate photons of frequency u is given 
by 
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(4) 



, 47r i^p 1 

where p ^ c/ (Sfi) is the curvature radius of the last closed 
line near the light cylinder in a dipole field. For v = 10^^ 
Hz, corresponding to 4 keV, we have: 



7~(-:/P)* = 8x10*. 
9 



(5) 



The radiative life time of the electrons is long enough to 
allow them to reach the light cylinder, namely: 



" I dj/dt I Stt Te v' 



(6) 



where is the classical radius of the electron; for P = 2.3 
ms and v = 10^^ Hz we have Tc — 5 s. These electrons 
could be accelerated either in the outer gaps or above the 
polar caps. In the latter case, we expect that a copious 
number of electron-positron pairs will reach, after suffer- 
ing radiation losses, the light cylinder with energy high 
enough for X-ray emission. 

An origin of the high-energy X-rays in a narrow outer 
gap can also explain the narrow profile (intrinsically < 
0.04 phase; Kuiper et al. 1998) measured in the ROSAT 
light curve, as well as the double peak profile, e.g. for 
the wide fan-beam geometry proposed by Romani (1996). 
However, an additional constraint is needed: the length of 



the gap over which the hard X-rays are produced should 

be small. In fact, the very strong curvature of the magnetic 
field lines would produce a broader pulse profile for a long 
gap. However, as Ho (1989) pointed out, strongly curved 
magnetic field lines, as in millisecond pulsars, enhance the 
production of non-thermal emission, allowing outer gaps 
shorter than those of normal radio pulsars. 

According Kuiper et al. (1999), PSR JG218+4232 is 
the most likely counterpart of the high-energy EGRET 
source 2EG J0220+4232 / 3EG J0222+4253 for energies 
between 100 MeV and 1 GeV. This would suggest a cut-off 
energy at least four orders of magnitude higher than the 
maximum energy of 10 keV to which the BeppoSAX data 
showed evidence for the detection of this source. Further 
high sensitivity observations at hard X-ray energies up to 
the high-energy 7-rays are therefore particularly relevant 
for the understanding of the physics of this and other mil- 
lisecond pulsars. 
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Fig. 6. PSR J0218-I-4232 spectra as detected by the MEGS. 
The lines represent the power-law used as fitting model 
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